Thermodynamics of a Trapped Bose-Fermi Mixture 
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By using the Hartree-Fock-Bogoliubov equations within the Popov approximation, we investigate 
the thermodynamic properties of a dilute binary Bose-Fermi mixture confined in an isotropic har- 
monic trap. For mixtures with an attractive Bose-Fermi interaction we find a sizable enhancement of 
the condensate fraction and of the critical temperature of Bose-Einstein condensation with respect 
to the predictions for a pure interacting Bose gas. Conversely, the influence of the repulsive Bose- 
Fermi interaction is less pronounced. The possible relevance of our results in current experiments 
on trapped *^Rb— ''"K mixtures is discussed. 

PACS numbers: PACS numbers:03.75.Kk, 03.75.Ss, 67.40.Db, 67.60.-g 



The recent experimental realization of ultracold 
trapped Bose-Fermi (BF) mixtures of alkali-metal atoms 
introduces an interesting new instance of a quantum 
many-body system 0, 12, ^-^i^ ^l^o stimulated a 
number of theoretical investieations that address, for ex- 
ample, the static property the phase diagram and 
phase separationja, stability conditions ^ and collec- 
tive excitations US 111 of trapped BF mixtures. These 
investigations have mainly concentrated at zero temper- 
ature using the standard Gross-Pitaevskii (GP) equation 
for the Bose gas, in which all the bosonic atoms are as- 
sumed to be in the Bose condensate. An extension of 
these theories at finite temperatures where the conden- 
sate is strongly depleted is therefore of high interest, and 
will also have practical applications. In the theory for 
a pure Bose gas, the simplest generalization of the GP 
equation including the effect of the noncondensed atoms 
in a self-consistent manner is the Popov version of the 
Hartree-Fock-Bogoliubov (HFB) approximation jl^. As 
discussed in Ref. |0, this approximation is expected to 
be good for both low and high temperatures. 

In this paper, we generalize the HFB-Popov approx- 
imation to binary BF mixtures and address the ques- 
tion of how the BF interaction affects the thermodynamic 
properties of mixtures. We calculate self- consistently the 
temperature-dependent density profiles of mixtures, as 
well as the condensate fraction and the critical tempera- 
ture of Bose-Einstein condensation (BEG), at various BF 
interaction strengths. Our present results provide the 
first self-consistent calculation of these thermodynamic 
quantities within the HFB-Popov theory which goes be- 
yond the semiclassical approximation used previously for 
determining the critical temperature IllL Il2| and den- 
sity profiles of binary BF mixtures [iJl- Our calcula- 
tions also show a highly nonlinear dependence of these 
quantities on the BF interaction. In the presence of the 
BF attraction, the thermal depletion of the condensate 
is remarkably decreased and the critical temperature is 
shifted towards high temperatures. Gonversely, the re- 
pulsive BF interaction affects the condensate fraction and 
critical temperature in the opposite direction. However, 



its influence is less pronounced compared to the attrac- 
tive case. 

The HFB-Popov mean-field theory for an inhomoge- 
neous interacting Bose gas has been derived in detail by 
Griffin in Ref. 10] . The generalization of this theory to 
a dilute binary BF mixture is straightforward. Here we 
merely give a brief summary of basic equations and em- 
phasize the necessary modification in the presence of the 
BF interaction. The trapped dilute mixture is portrayed 
as a thermodynamic equilibrium system under the grand 
canonical ensemble whose thermodynamic variables are 
Nb and Nf, respectively, the total number of trapped 
bosonic and fermionic atoms, T, the absolute tempera- 
ture, and /if, and /iy, the chemical potentials. The density 
Hamiltonian of the system is given by (in units of h = 1) 



= Tib + Ti-f + Ti-bf, 

,2 

Hb - yj+{r) ' 



2nib 



2m 



f 



Hf = 0+(r) 

nbf = .g6/V+(r)^(r)0+(r)0(r), 
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where ip{r) (0(r)) is the Bose (Fermi) field operator 
that annihilates an atom at position r. Here we con- 
sider a spherically symmetric system, with trap poten- 
tials V^t'ap('') = ^^bji^b f'''^ where rubj are the atomic 
masses and cubj are the trap frequencies. The interac- 
tion between bosons and between bosons and fermions 
are described by the contact potentials and are param- 
eterized by the coupling constants gbb = 4:nh'^ abb / nib 
and Qbf — inh^abf /rur in terms of the s-wave scattering 
lengths abb and abf, with = nibm f / {rrib+m f) being 
the reduced mass. We neglect here the fermion-fermion 
interactions, since we are considering a spin-polarized 
Fermi gas where s-wave collisions are forbidden by the 
Pauli principle. In the dilute regime, we may treat the 
density Hamiltonian describing the BF coupling in a self- 
consistent mean-field manner, namely, 

Hbf ^ 9bf [^-+^-(0+0) + (V+V'>0+'^- • 
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This kind of decomposition has been used extensively for 
theoretical investigations of BF mixtures at zero temper- 
ature p^ . 

To the Bosc field operator ip{r,t), we shall apply 
the usual decomposition into a c-number part plus an 
operator with vanishing expectation value: 'ip{r,t) = 
<i>(r)e~*'^''~^''''* + ^{r) represents the condensate 

wave function with eigenvalue eg and the operator ip{r) 
represents the excitations of the condensate. This ansatz 
is then inserted in the equation of motion for iplr^t): 



.dip 



V 

2mb 



* trap 



(2) 

The statistical average over Eq. ((2Jl and the replacement 
of the cubic term ^^iptp by the average in the mean-field 
approximation 2{'ip'^ip)'ip with neglecting the anomalous 
expectation value {'ip'<P) s-nd its complex conjugate lead 
to the generalized GP equation, 



CgpHv) = 0, 



(3) 



where £gp = -\/'^/2mb + T^^-ap - + 9bb{nc{r) + 
2n{r))+gbfnf{r) with the local density of the condensate 
nc(r) = |<i>(r)|^, of the depletion ri(r) — (-0+(r, t)?/'(r, t)), 
and of the Fermi gas ri/(r) — ((/i+(r, i)0(r, t)). The 
condensate wave function in Eq. (O is normalized to 
= l/{ef^(^o-i^'') - 1) with 13 = {ksT)-^. 
The subtraction of Eq. J^J from Eq. (O gives rise to 
two coupled equations of motion for iplr^t) and its ad- 
joint, which can be solved by the usual Bogoliubov trans- 
formation, ^{v,t) = [ui(r)die~"'* -I- w*(r)(3;^e"'*] , 
to the new Bose operators and af . This gives the 
coupled Bogoliubov-dcGcnnes (BdG) equations, 

[Cgp + gbbTicir)] Ui(r) + gbbnc{Y)vi{Y) = eiUi(r), 
[Cgp + gbbnc{Y)]vi{v) + gbbnc{v:)ui{v) = -ejWi(r).(4) 

These equations define the quasiparticle excitation ener- 
gies Cj relative to the condensate eigenvalue £o, and the 
quasiparticle amplitudes Ui and Vi. Once these quantities 
have been determined, the density of the depletion is ob- 
tained in terms of the thermal number of quasiparticles 
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V27r2 \ zel3'-\ 2 2eHF J 

X {euF - V^ttap + £o - '2-gbbnb{v) - gbfnf{v)Y'^, 



where z = e^^^-^-^"^^ = 1 + = {e^ + glbnl{Y)Y''^ 

and nb(r) — nc{v) + n(r) is the total density of the Bose 



gas. In the above equations, to eliminate the numerical 
errors due to the necessary truncation of the numerical 
basis set, we adopt the strategy of Ref. 0| and introduce 
an energy cutoff E^, above which the semiclassical local- 
density approximation has been employed. 

To solve the generalized GP and BdG equations, one 
has to find the local density of the Fermi gas ?^/ (r) . To 
this end, we insert 0(r, t) = (pi{r)cie~'^'^'* in Eq. to 
diagonalize the quadratic Hamiltonian for (/)(r, t) in terms 
of the new Fermi operator Ci that annihilates a fermion 
at state tpi{r). This leads to a stationary Schrodingcr 
equation for (^i(r). 



+ Wrap+ gbfin-c + n) 



2m/ 

The density of the Fermi gas is thus obtained by 
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where {c^Ci) — (e'^^''^^-''^ -f 1)^-'^ is the Fermi distribu- 
tion. Analogously to Eq. (O we have applied the finite- 
temperature Thomas-Fermi (TF) approximation only for 
high-lying Fermi levels above an energy cutoff El to 
avoid the truncation errors. 

Equations ©-fZl form a closed system of equations 
that we have referred to as the "HFB-Popov" equations 
for a dilute BF mixture. We have numerically solved 
these equations by an iterative procedure as follows: The 
generalized GP and BdG equations are first solved self- 
consistently for $(r), Ui(r), and Wi(r) as described in Rcf. 
[Tol l to evaluate nc(r) and n(r), with n/(r) set to the re- 
sult for an ideal Fermi gas. Once nc(r) and n(r) are 
known, the eigenfunctions in Eq. @ are obtained nu- 
merically and are used to update n-/(r) in Eq. ((TJ. This 
newly generated «/(r) is then inserted in the GP and 
BdG equations and the process is iterated to convergence. 
At each step, the chemical potentials for the Bose gas and 
the Fermi gas are fixed by the normalization conditions, 
J c?rnfc(r) = Nb and J drnf{r) = Nf, respectively. 

As an illustration of this procedure, we consider a 
mixture of 2000 ^'^Rb (boson) and 2000 4°K (fermion) 
atoms in an isotropic harmonic trap, for which the order 
parameter $(r), the quasiparticle amplitudes Ui{r) and 
Vi{r), and the orbits (pi{r) can be labelled by {n,l,m), 
according to the number of nodes in the radial solution 
n, the orbital angular momentum I, and its projection 
TO. In addition, we use the following parameters j^: 
TOfc = 1.45xl0~25 kg, LJb = 27rx216 Hz, to//to6 = 0.463, 
ujf/ujb = 1.47, abb — 99ao, and abf = — 410ao, where 
flo = 0.529 A is the Bohr radius. Because our calculations 
are especially delicate near the critical temperature, we 
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(a) T = 80nK (b) T = 110 nK 
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FIG. 1: The density profiles of the condensate nc(r) and of 
the Fermi gas nF(r) for a mixture of 2000 '^'^Rb and 2000 ""K 
atoms in an isotropic harmonic trap at various BF interaction: 
a;,/ =0 (full line), atf = +410ao (dashed line), and atf — 
— 410ao (dash-dotted line), where ao = 0.529 A is the Bohr 
radius. Insets show the density profiles of the noncondensate 
n(r). We have taken atb = 99ao and = 2n x 216 Hz. 
The coordinate r and densities are measured in units of the 
harmonic oscillator length a^^ and (a^^^ , respectively. 



have taken rimax = 32, /max = 64 and high energy cutoffs 
of = GOhoJb and E[ = QOfiujb to ensure the accuracy. 
Throughout the paper, we also express the lengths and 
energies in terms of the characteristic oscillator length 
= {'h/miiijJbY^'^ and characteristic trap energy fiLOb, 
respectively. 

In Fig. 1, we present our results for the density pro- 
files of the condensate, of the noncondensate, and of the 
Fermi gas at two temperatures. The cases with the BF 
interaction and without the BF interaction are shown by 
the dash-dotted lines and full lines, respectively. We have 
also considered a fictitious case of a positive BF interac- 
tion: a-bf = +410ao (dashed lines). The choice of the 
first temperature, T = 80 nK, corresponds to the situa- 
tion in which the condensate and noncondensate have an 
approximately equal number of atoms, while the other 
temperature T = 110 nK is chosen to be close to the 
critical temperature for a pure interacting Bose gas with 
the same number of bosons, Tc « 112 nK. As clearly 
emerges from the figure, the density profiles of the con- 
densate and of the Fermi gas are strongly affected by the 
BF interaction at both temperatures. In particular, the 
densities around the center are significantly enhanced in 
the case of the BF attraction. The density profile of the 
noncondensate (shown in the insets), on the other hand, 
is less influenced by the BF interaction due to its broad 
distribution and the strong repulsion from the conden- 
sate. 

In Fig. 2, we show our predictions for the temperature 
dependence of the condensate fraction Nc{T)/Nb. The 
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FIG. 2: Temperature dependence of the condensate frac- 
tion. The dashed and dash-dotted lines correspond to atf = 
-|-410ao and — 410ao, respectively. The solid lines are the 
result for a pure interacting Bose gas. The empty cir- 
cles are ideal Bose gas result with the finite size correc- 
tion, N°/Nb = 1 - {T/T'^f - 2.1825(^/^°)27V~^''^ where 
ksT^ ~ 0.94hujtN^^'* . The inset highlights the condensate 
fraction near the critical temperature. The other parameters 
are the same as in Fig.l. 



essential feature of the figure is the importance of the at- 
tractive BF interaction that results in a sizable quenching 
of the thermal depletion compared to the prediction for 
a pure interacting Bose gas. Contrarily, the effects of the 
repulsive BF interaction are more subtle and are always 
very small. The sizable enhancement of the condensate 
fraction predicted by our calculation follows from the fact 
that in the presence of the BF attraction the conden- 
sate effectively experiences a more tightly confining po- 
tential. As a consequence, if we neglect the corrections 
due to the interaction between bosons and the finite size 
effect, the critical temperature T° = 0.9AfiujeffNl^^/kB 
is effectively increased and the condensate fraction is, 
therefore, enhanced according to the ideal gas result 
NjNb = 1 - (r/T0)3. 

Closely related to the condensate fraction, another im- 
portant parameter characterizing the effect of the BF in- 
teraction is the shift of the critical temperature from the 
pure interacting Bose gas case. In Fig. 3, we report the 
HFB-Popov results for the relative shift of the critical 
temperature ST^/Tc as a function of a^/ in solid circles. 
Here is determined as the maximum of the function 
(PNc/dT'^ nil. The semiclassical predictions for ST^/T^, 
calculated as in Ref. Tl'l in the first order of are 
also shown by the dashed line. The agreement of these 
two approaches is reasonably good for a weak BF inter- 
action {\abf\ ^ lOOao). However, as \abf\ increases, our 
HFB-Popov results diverge from the semiclassical predic- 
tions. In particular, for the realistic BF s-wave scattering 
length for *^Rb— '^''K mixtures, abf = — 410ao, the devi- 
ation becomes remarkable. 

We now turn to consider the experimental relevance of 
our results. In current experiments, the realistic number 
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FIG. 3: The relative shift of the critical temperature as a 
function of the BF s-wave scattering length ai,f (in units of 
ao). The solid circles are the result calculated using the HFB- 
Popov equations and the full line is a parabolic fit to guide eye. 
The dashed line is the result calculated by Eq. (18) in Ref. 
[11]. The inset shows the functions dNc/dT and SNc/dT"^ for 
the case of a;,/ = — 410ao. Tc is extracted from the maximum 
of Nc/dT"^ . The other parameters are the same as in Fig.l. 




a few TiLOb and applying the TF approximation for the 
whole Fermi spectra. The accuracy of this semiclassical 
treatment has been checked by the comparison with the 
full quantum mechanical calculations for a small mixture. 
The condensate fraction obtained by these two methods 
coincides within 1% errors. 

In Fig. 4, we present the results for Nc{T)/Nb of a 
mixture of 2 x 10** ^''Rb and 5 x 10* '^''K atoms con- 
fined in an isotropic trap with Wf, = 2ti x 91.7 Hz. This 
choice corresponds to the typical experimental situations 
at LENS |23|. Although the trap becomes more shal- 
low, the enhancement of the condensate fraction shown 
in Fig. 4 is quantitatively similar to the one of Fig. 2, 
due to the much larger values of and Nf contained 
in the trap. Finally, in this case we also roughly esti- 
mate the relative shift of the critical temperature due 
to the BF attraction, {6Tc/Tc)bf — +4%, which is com- 
parable to the shift due to the boson-boson interaction. 



-3.2%, and the fi- 



-1/3 _ 



-2.7% 



nite size correction, {STc/Tc)fs — ~'0■73N^ 

M 

In conclusion, we have generalized the HFB-Popov 
theory to binary BF mixtures and have presented a de- 
tailed study of the thermodynamic properties of mixtures 
at finite temperature, including the density profiles, the 
condensate fraction, as well as the critical temperature 
of BEC. These quantities are found to depend on the 
BF interaction in a nonlinear way. Moreover, under 
conditions appropriate to the ^^Rb— ''^"K mixture in the 
LENS experiments, the condensate fraction and the 
critical temperature of BEC are significantly enhanced 
with respect to the prediction for a pure interacting Bose 
gas. This enhancement might be observable in current 
experiments. 



FIG. 4: Temperature dependence of the condensate fraction 
for a mixture of 2 x 10'' *^Rb and 5 x lO* ""^K atoms confined 
in an isotropic trap with uui, = 2Tr x 91.7 Hz, calculated in the 
semiclassical version of the HFB-Popov theory as mentioned 
in the text. 



of ^^Rb and ^''K atoms in mixtures is about ten times 
larger than what we assumed here Q. For such large 
number of atoms, our calculation is very time consum- 
ing and we then resort to the semiclassical version of 
the HFB-Popov theory by setting to an energy of 
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